Astronomy & Astrophysics manuscript no. stellar'noise'and'planetary 'mass'detection'limits © ESO 2010 

October 14, 2010 



Planetary detection limits taking into account stellar noise ^ 
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ABSTRACT 



Context. Stellar noise produced by oscillations, granulation phenomena (granulation, mesogranulation and supergranulation) and ac- 
tivity affects radial velocity measurements. The signature of this noise in radial velocity is small, around the meter-per-second, but 
already too much for the detection of Earth mass planets in habitable zones. 

Aims. In this paper, we address the important role played by observational strategies in averaging out the radial velocity signature of 
stellar noise. We also derive the planetary mass detection limits expected in presence of stellar noise. 

Methods. We start with HARPS asteroseismology measurements for 4 stars (fiHyi, o-CenA, yu Ara and rCeti) available in the ESO 
archive plus very precise measurements of aCenB. This sample covers different spectral types, from G2 to Kl and different evolu- 
tionary stage, from subgiant to dwarf stars. Since the span of our data ranges between 5 to 8 days, only stellar noise sources with 
a time scale smaller than this time span will be extracted from these observations. Therefore, we will have access to oscillation 
modes and granulation phenomena, without important contribution of activity noise which is present at larger time scales. For those 
5 stars, we generate synthetic radial velocity measurements after fitting corresponding models of stellar noise in Fourier space. These 
measurements allows us to study the radial velocity variation due to stellar noise for different observational strategies as well as the 
corresponding planetary mass detection limits. 

Results. Applying 3 measurements per night of 10 minutes exposure each, 2 hours apart, seems to average out most efficiently the 
stellar noise considered. For quiet KIV stars as o-CenB, such a strategy allows us to detect planets of ~ 3 times the mass of Earth 
with an orbital period of 200 days, corresponding to the habitable zone of the star. Our simulations moreover suggest that planets 
smaller than typically 5 can be detected with HARPS over a wide range of separations around most non-active solar type dwarfs. 
Since activity is not yet included in our simulation, these detection limits correspond to a case, which exist, where the host star has 
few magnetic features. In this case stellar noise is dominated by oscillation modes and granulation phenomena. For our star sample, a 
trend between spectral type and surface gravity and the level of radial velocity variation is also emphasized by our simulations. 

Key words, stars: individual: /JHyi - stars: individual: yu Ara~ stars: individual: o-Cen A - stars: individual: rCeti - stars: individual: 
cf CenB - stars: planetary systems - stars: oscilation - techniques: radial velocities 



1. Introduction 

The majority of planets discovered by Doppler spectroscopy 
are gaseous giants similar to Jupiter, with masses of a hun- 
dred times the mass of the Eartl(3(Me). However, since approx- 
imat ely 5 years, planets fr om 2 to 10 Mm, h ave been detected 
(e.g. iMavor et al.ll2009blli lUdrv et al.ll2007l) . This has become 
possible thanks to the stability and precision of new generation 
high-resolution spectrographs, and because of dedicated obser- 
vational strategies allowing to average out p erturbations coming 
from stellar oscillations (Sa ntos et 



age out p e 
al] l2004 . 



Pressure waves (p-modes) propagate at the surface of solar- 
type stars leading to a dilatation and contraction of external en- 
velopes over time scales of a few minute s (5 to 15 minut es for 
the Sun; Sch riiver & Zwann 2000; Broom hall etani2009l) . The 
radial-velocity signature of these modes is typically varying be- 
tween 10 and 4( 30 cm s~', depending on t he star type and evo- 
lutionary stage dSchriiver & Zwannll2000l) . The amplitude and 



period of oscillation modes increases with mass along the main 
sequence. Theory predicts that frequencies of p-modes rise up 
with the square root of the star mean density and that their 
amplitudes are proportional to the luminosity over mass ratio 
(Christensen-Dalsgaard 2004). The most precise spectrograph 
nowadays, HARPS, can reach a radial-velocity precision better 
than the meter-per-second in a short exposure time for bright 
stars (typically 1 m s~' in 1 minute for a V = 7.5 K dwarf, 
iPepe etal.l2005 h. These acoustic modes are then directly observ- 
able and can mask an eventual smallest planet signature. 

The different phenomena of granulatiorQ (granulation, meso- 
granulation, and supergranulation), due to the convective nature 
of solar type stars, also affect radial velocity (RV) measurements. 
These convective phenomena can be found over all the stellar 
surface, exc ept in active regions where convection is greatly re- 
duced (^e.g iDravinsI 1 1 9821; Livingstonl 1 1 9821: lBrandt & Solankil 
ll990l;lGravlll992 : iMeunier et al .1120101) . The ampUtudes of gran- 



* Based on observations collected at the La Silla Parana Observatory, 
ESO (Chile), with the HARPS spectrograph at the 3.6-m telescope. 
' see The Extrasolar Planets Encyclopaedia, http://exoplanet.eUj 



- hereafter, the term granulation phenomena will be used to describe 
all type of convective motion (granulation, mesogranulation and super- 
granulation), whereas granulation will just be used for the smaller time 
scale convective motion. 
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ulation phenom ena are similar to the ones observed for pres - 
sure modes (e.g.'Schriiver & Zwanij200dllKjeldsen et al.l2005l) . 
Granulation corresponds to a small convective pattern with a 
lifetime shorter than 25 minutes, and a diameter smaller than 
2 Mm (.Title et al.lll989llDel Moroll200l . At much larger scale, 
we can find supergranulation. This phenomenon, linked to very 
large convective patterns, 15 to 40 Mm, can have a lifetime up 
to 33 hours in the Sun (Del Moro et al. 2004). Mesogranulation 
represent a convective phenomenon which can be located be- 
tween the granulation and supergranulation, in terms of size 
and lifetime (Harvev 1984; Palle et al. 1995; Sc hriiver & Zwaniil 
12000 ). At the moment, the HARPS-GTO strategy, for the very 
high precision programs, uses long exposure times (15 minutes), 
in order to reduce the effects produced by stellar oscillation on 
RV measurements. However, this does not reduce low frequency 
noises coming from granulation phenomena. 

On the time scale of tenths of days, another type of noise, 
due to the presence of activity related spots and plages, may 
perturb precise RV measurements. The presence of spots and 
plages on the surface of the star will break the flux balance be- 
tween the red-shifted and the blue-shifted halves of the star As 
the star rotates, a spot, or a plage, move s across the steUar disk 
and produces an apparent Doppler shift (ISaar & Donahudl 19971; 
lOueloz et al.ll2001l; iHuelamo et alj|200"8t iLagrange et alj|2010h . 
This modulation can be hard to disentangle from the Doppler 
modulation caused by the gravitational pull of a planet. On the 
Sun, taken as a proxi for G stars, the amplitude of this pertur- 
bation can reach 40 cm s ' in high activity phase dMeunier et alJ 
l20Toh . Spots are cooler than the solar mean surface temperature 
and plages hotter Since active regions contain both, the noise 
induced by spots and plages will usually compensate, but not 
entir ely since the surface ratio between s pots and plages varie s 
(e.g. lChapman et alj|200lh . According to Meunier et al] (1201 Oh . 
the major effect of activity is not due to the break of the flux bal- 
ance but to the inhibition of convection in active regions. They 
found that the corresponding noise could vary between 40 cm s"' 
at minimum activity and 140 cm s"' at maximum. 

In this paper, we present simulations exploring new measure- 
ment strategies allowing to average granulation phenomena and 
oscillation modes at the same time. The detection limits for dif- 
ferent strategies are derived, pointing out in particular a strategy 
that efficiently optimizes detection ability and realistic observa- 
tional cost. Our simulations also emphasize that for dwarf stars, 
the level of RV variation, mainly due to stellar noise, is corre- 
lated to the spectral type and/or the surface gravity, log g. 

We note that in this study we will just focus on understanding 
the stellar noise influencing in time scales smaller than the span 
of our data, 8 days. On such time scales, only the influence of os- 
cillation and granulation phenomena can be fully characterized. 
Longer term activity signals, typically modulated at timescales 
of ~30days (the usual rotation period of solar type stars) are 
not fully accessible. However, noise coming from activity re- 
lated phenomena are still present in our data. The obtained re- 
sults will thus only be valid for stars without significant activity 
phenomena. The effect of plages and spots will be estimated in 
a forthcoming paper 

2. Data selection 

In order to characterize as well as possible the RV noise in- 
trinsic to Sun-like stars, we wanted the most precise asteroseis- 
mology measurements (continuous follow up of the star with a 
high sampling rate). Another requirement was to use only data 
from one instrument, to avoid any instrumental discrepancies. 



We thus choose to use only HARPS measurements, since this 
spectrograph can reach a RV precision better than the meter per 
second, on the long term (Pepe et al. 2005), and even a better 
precision, at the lev el of a few 10s of cms ', over one night 
dBouchv et al.ll2005l) . Four stars with adequate observations are 
available in the HARPS ESO archive: ySHyi (HD2151), nAm 
(HD 160691), a Cen A (HD 128620) and t Ceti (HD 10700). The 
measurements of aCenB (HD 128621) have been provided by 
Pepe, F. (private communication) and do not correspond to con- 
tinuous measurements. This star has been observed for 8 con- 
secutive nights, though only 3 blocks of 20 minutes (about 27 
measurements of 13 seconds) were done during each night. We 
will see afterwards that this is sufficient to characterize all type 
of considered noises. Depending on the star, the total span of the 
measurements varies from 5 to 8 nights. Table[T]gives more in- 
formations about the stellar parameters. Our sample contains 3 
main sequence (MS) stars, 1 sub-gia nt and I s tar, fi Ara, which 
is between the 2 regimes ( S oriano & Vauclaiij|201Q) . This will 
allows us to study a possible correlation between the level of RV 
variation and spectral type and/or evolution. 

Exposure times are always shorter than 100 s. Due to CCD 
readout, a dead time of ~30 s exists between each exposure. Such 
a strategy defines precisely the structure of oscillation modes, 
since they have a typical time scale above 5 minutes. The con- 
tinuous follow up of the star during 5 to 8 days is essential in 
order to characterize convective motions (granulation phenom- 
ena), which have time scales up to ~30 hours. In the case of 
a Cen B, the 20 minutes of continuous measurements allow us to 
characterize granulation, since it has typical time scales smaller 
than 25 minutes. The mesogranulation perturbations can be ana- 
lyzed thanks to the repetition, 3 times a night, of the continuous 
measurements. And finally, the 8 consecutive nights allow us to 
address supergranulation. We will see in Sec. 13.21 that for the 
case of a Cen B, all types of noise can be distinguish clearly in 
the power spectrum, as it is the case for the other stars with con- 
tinuous measurements. 

In our sample, // Ara is the only star known to have plan- 
etary companions. We clearly see in the RVs the effect of the 
small planet in the system, which we suppress taking the orbital 
parameters given in lPepe et"al1 (l2006l) . 



3. Generating synthetic RV measurements 

In this paper, we study different observational strategies to av- 
erage out as much as possible the effect of stellar noise on RV 
measurements. The main goal is to find a strategy capable of de- 
tecting very small mass planets in habitable regions. To be able 
to search for such long period planets, observational strategies 
must be efficient on several years. Since the asteroseismology 
data of the selected stars span only over a maximum of 8 days, 
we have to create longer synthetic RV data, containing the same 
structure and level of noise. To generate these synthetic data, we 
first calculate the velocity power spectrum density (VPSD) us- 
ing the periodogram technique. We then fit the obtained spec- 
trum with a noise model and go back to the RV space using 
ffie property that the periodogram is eq uival ent to fitting series 
of sine waves (e.g. Scargle 1982HZechmeister & Kursteili2009l) . 
The transformation from RVs to the VPSD, as well as the in- 
verse one, does not loose any information, so the synthetic RV 
sets should present the same noise structure as the real data we 
started from. We have however to keep in mind that even if we 
have synthetic RV data over more days then the real measure- 
ments, we do not take into account, in our model, perturbing 
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Table 1. Properties of ot ir sample stars. Fo r ea ch star, the dead tirn e d ue to CCD readout is 31s ((a) ISousa et al ] (120081) . (b) 
iBedding et all (l2007h . (c) iBazot et all (l2007h . rdl lBouchv etaP (l2005h . (e^ lTeixeira etaP (l2009h V ST is for spectral type and T„p 
for exposure time. 



Star 


Mv 


Mass [Mj"' 


T,ff iKr> 


[Fe/Hf> 


log 


log(R'«jf) 


ST 


# meas. 


# nights/span [days] 


T„p [s] 


/?Hyi 


2.8 


1.10 


5872 


-0.11 


3.95 




G21V 


2766'*' 


6/6 


40 to 50"" 


yu Ara 


5.2 


1.08 


5780 


0.30 


4.27 


-5.20 


G3IV-V 


2104('" 


8/8 


100*''* 


Q-Cen A 


0.0 


1.09 


5623 


0.22 


4.22 


-5.07 


G2V 


4959(f) 


5/5 


2 to lO*'^' 


rCeti 


3.5 


0.63 


5310 


-0.52 


4.44 


-4.93 


G8V 


1962*''' 


5/6 


23 to 40<"' 


QfCenB 


1.3 


0.91 


5248 


0.26 


4.55 


-4.98 


KIV 


783 


8/8 
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effects with time scales longer than the real observation time. 
Therefore activity noise is not fully included. 

3. 1. Velocity power spectrum density 

The first step is to calculate the velocity power spectrum of the 
asteroseismology RV measurements. Since the measurements 
are not evenly spaced in time, we can not use a standard FFT 
algorithm. We thus use a weighted Lomb-Scargle periodogram 
(^mb 1976; Scargle 1982; Zechmeister & Kurster 2009) to cal- 
culate the power spectrum of the HARPS RV measurements. The 
minimum frequency of the spectrum is set to 1/T, where T is the 
total span of the observing run and the maximum frequency is 
set to the Nyquist frequency for the median sampling. To obtain 
the VPSD, which is independent of the observing window, we 
multiply the power by the effective length of the observing run, 
calculated as the inver se of the are a under the spectral window 
(in power, see Kjeldse net al.ll2005l) . 

In the raw VPSD of all the stars, except cCenB, we no- 
tice huge excess of power at 3.1, 6.2 and 9.3 mHz (see Fig. [T] 
top panel). It as been shown that these anomalies are due to a 
problem in the guiding engine of the telescope, resulti ng in a 
poor guiding prec ision (e.g Teixeira et al. 2009; Beddin g et alJ 
120071; iBazotet alj|2007l: Carrier & Eggenberger 2006). The data 
for cCenB, taken in 2009, does not present these peaks due to 
corrections in the guiding software. 

To correct this noise effect in our measurements, we use 
the method devel oped for a Ce n A ( B utler et al. 2004), ySHyi 
(iBedding et al. 2007) and rCeti (iTeixeira et al. 2009). The idea 
is to down weight bad points produced by the poor guiding. First 
of all, we compute the raw power spectrum and remove the high- 
est peaks corresponding to oscillations or granulation phenom- 
ena. We subtract the related sinusoid component from the RV 
data and we restart the process of calculating the periodogram 
and selecting the highest peak until all the power coming from 
real signals have been killecQ. The residual velocities obtained 
should only be composed of noise. If this noise is pure photon 
noise, the residual velocities should respect a gaussian distribu- 
tion. For each night, we select all the point deviating from the 
mean scatter by more than 3cr and we multiply their statistical 
weight by the probability of the point to belong to a gaussian 
distribution. Finally, we scale all the statistical weights, night by 
night, to reflect the noise measured at high frequencies. Figure 
[1] shows, for a Cen A, the improvement brought by the cleaning 
process. After down weighting bad points, the noise level at high 
frequencies is reduced by 25 to 30 % for all stars, except a Cen B 
that do not present this guiding noise effect in the raw data. In 



^ For yu Ara, a Cen A and cfCenB, all the peaks belo w 2.5, 3.5 and 
6 mH z are removed, respectively. For rCeti, following ITeixeira et alJ 
( l2009h , we kill all the frequencies below 0.8 mHz and between 3.1 and 
5.5 mHz. 
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0.006 0.008 
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0.010 0.012 



Fig. 1. Cleaning of the guiding noise for ffCenA. Top 
panel. Powei spectrum of the RV measurements, using the raw 
errors to calculate the statistical weights. We can see clearly the 
peaks at 3.1, 6.2 and 9.3 mHz due to the guiding noise. Bottom 
panel.Same power spectrum but after the cleaning process. The 
guiding noise peaks are much fainter 



the case of /3 Hyi, the power spectrum has been kindly provided 
by Tim Bedding and without having access to the raw RVs, we 
could not correct for this guiding noise effect (the peaks at 3.1 
and 6.2 mHz are visible on the jSHyi power spectrum in Fig.|2]i. 
The corrected VPSD will be used for the following of the study. 

For all the stars in our sample, we observe a raise of the 
VPSD towards low frequencies (see Fig.|2]) due to the nature of 
granulation phenomena. The decrease with frenquency follows 
until reaching the Lorentzian p-modes bump. Finally, after the 
signature of oscillation modes, we can see that the amplitude of 
the VPSD stays at the same level due to photon and instrumental 
noise. 



3.2. Fitting tlie veiocity power spectrum density 

Once the VPSD is calculated, we need to fit it to obtain a spec- 
trum that is independent of the original data sampling. This will 
also allow us to separate individually each type of noise, which 
will give us an idea about the RV contribution of each of them. 

The part of the VPSD coiTesponding to granulation phenom- 
ena can be explained by a model in which each source of con- 
vective motions i s desc ribed by an empir ical law initia l ly pro - 
posed bv lHarvevl(ll984l) and reviewed bv iAndersen et al.l (Il994l) : 
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Fig. 2. VPSDs and fits for the 5 stars of our samples. The grey dots corresponds to the binned spectrum using a boxcar algorithm 
with a changing width according to the frequency. The global fit including granulation phenomena, oscillation modes, photon noise 
and instrumental noise (black line) is obtained on the binned spectrum. Each contribution to the fit is shown. We can see the 3 
types of granulation fitted in light grey lines (supergranulation, mesogranulation and granulation from left to right respectively). 
The Lorentzian fit to p-modes corresponds to the dashed line and the constant, fitting the photon and instrumental noises, can be see 
in dotted fine. 



IPaile et al.l (Il995h . This model corresponds to an exponentially 
decaying function with frequency: 

where P(v) is the VPSD, A the power density of the correspond- 
ing convection motion, B its characteristic time scale and C the 
slope of the power law. The power remains approximatively con- 
stant on time scales larger than B and drops off for shorter time 
scales. The total VPSD due to granulation phenomena is the sum 
of the VPSD for each phenomenon (G = granulation, MG - 
mesogranulation and SG = supergranulaion): 

Pto,(v) = Pc(v) + Pmc(v) + Psciv). (2) 



The frequency of the oscillation modes (p-modes) can ana- 
lytically be derived using a perturbation theory and the asymp- 
totic theory dTassoulll 19801) . leading to: 

v„, ~ Av^n + ^ -t-ffj, (3) 

where Av is the large separation, n the radial order and / the an- 
gular degree of the mode. We are interested in the total power 
due to p-modes and not in e ach mode separately . We thus 
use the technique proposed by Kieldsen et al.l (12005 ^ (see also 
Arentoft et al. 2008) to produce a single hump of excess power 
that is insensitive to the fact that the oscillation spectrum has 
discrete peaks. This is done by convolving the observed VPSD 
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with a Gaussian having a full width at half maximurr|j of 4Av. 
Finally, we fit the convolved curve with a Lorentzian function 
(e.g. ,Lefebvre et al..,2008.) : 



P(v) = Al 



(v-V0)2+r2' 



(4) 



where Ai correspond to the amplitude of the Lorentzian, F to the 
full width at half maximum (FWHM) and vq to the mean of the 
Lorentzian. 

At very high frequency, after the p-modes bump, we can see 
the contribution of photon and instrumental noises. Supposing 
that these noises are gaussian, which is at least the case for 
photon noise, their contribution correspond to a constant power 
ranging over all the spectrum. 

To summarize, the entire VPSD can be adjusted by the fol- 
lowing function with 13 free parameters: 



/=1 



+ A, 



(v-vo)2 + F2 



■ constant. 



(5) 



where i correspond to each type of granulation phenomena. In 
order to perform the VPSD adjustment with 13 free parameters, 
we bin the spectrum, using a boxcar algorithnfl and use an itera- 
tive approach. We first adjust the Lorentzian and the constant on 
the oscillation modes, and then, fit the Harvey functions on the 
low frequency part of the spectrum, taking into account the con- 
stant derived before. Finally, we derive a global simultaneous fit 
of the model including the Harvey functions, the Lorentzian and 
the constant using as initial conditions the results of the individ- 
ual fits. The VPSD and final adjustment for each star can be seen 
in Fig.|2]and the different parameters, in Table|2] 

For the case of a Cen B, the observation data we have does 
not correspond to continuous follow up of the star. To check if 
it does not affect high frequency noises, we compare the high 
frequency spectrum of our data with one night of asteroseismol- 
ogy done on aCenB in 2004 (HARPS ESO archive). We find 
the same p-modes envelope in both case, which proves that the 
data we have for a Cen B are as good as the one gathered for the 
4 other stars (see Fig.O. 

Concerning oscillation modes, all the stars of our sample 
have been previously studied in details. We thus can compare 
the mean frequency of the p-modes bump fitted in this paper, 
with the one found in previous studies ( Arentoft et al. 2008; 
iKiel dsen et al. 2008, and reference therein). According to Fig. 
1 1 in .Arentoft et aL (.2008i) . the mean frequency of the p-modes 
bump are the following: 1 mHz for jSHyi, 2mHz for /vAra, 
2.5 mHz for a Cen A, 4.5 mHz for r Ceti and 4 mHz for a Cen B. 
These frequencies can be compared with our results, see vq val- 
ues in Table|2l and we find a very good agreement. 

To compare the fitted parameters of granulation phenom- 
ena with the literature values is more complex. To our knowl- 
edge, only t he Sun has been prev iously studied in detail 
(ILefebvre et a l. 2008 : Paile et alJll995i) . In addition, the normal- 
ization of the algorithms used by these authors to calculate the 



The value of Av for e ach star of our sampl e is taken from the liter- 
ature:57.2 /iHz for ySHyi jBedding et alll2007l) . 106 fiHz f or o-CenA 
dKieldsen et al. 200 5]), 90 nH z for u Ara dBazot et al.ll2005h. 169 uHz 
for T C eti ( fXeixeira et alj2009l) and 162 fiHz for a Cen B ( IKieldsen et alj 
12003) 

The width of the box changes according to the frequency. The goal 
is to obtain a binned spectrum which has nearly the same spacing be- 
tween the points in a logarithmic scale, which makes fit convergence 
faster. 



VPSD is different from the one adopted here. This implies the 
derivation of different amplitudes for the considered granula- 
tion phenomena. Some work could be done to normalize the 
values, but this is out of the scope of the present paper We 
can nevertheless com pare the different time scales obtained. In 
ILefebvre etal] (|2008|) and lPalle et~al] (Il995h . we have access the 
the characteristic time t, which is by definition equals to B/(2n). 
Changing the characteristic times from t to B, Lefebvre et al.l 
( l2008h obtain Bg = 22 min whereas Pall e et all ( Il995h get Be 
= 39 min. In our case, for a Cen A (G2V) and ju Ara (G3IV-V) 
which are close to the Sun in terms of spectral class, we find Be 
- 18 min and 44 min respectively, which is in good agreement. 
For the case of meso- and sup ergranulation, we find some dis- 
crepancies in the literature. In iPalle et alJ (fl99l . Bmg ~ nh 
and Bsc ~ 175 h. According to more precise an d recent results, 
we can extract roughly from Fig. 2 in lLefebvre~ et al. (2008) that 
Bsc ~ 2.10'^ Hz ~ 14/1. Notice that our results for aCenA 
(G2V) and /u Ara (G3IV-V) are close to this second value. 



3.3. Synthetic radial velocity measurements 

Once the VPSD fit is obtained, we have to calculate the RV for 
each date of an observational calendar Since the periodogram 
techn i que is equivalent to a w eighted sine wave fitting (Scargl^ 
119821: IZechmeister & Ktirsteiil2009b . to pass from the spectrum 
to RVs, we have to calculate for the desired epochs f,: 



RV{ti) = VVPSD(v)(sin(27rvf; + phase{v))). 



(6) 



where the frequency v goes from 1/T to the Nyquist frequency 
with a step of 1/T. Supposing that all the noises due to stellar 
perturbations are independent, we choose a uniformly random 
phas^ between and 2;r. 

To check if the computation of synthetic RV measurement 
over long period does not introduce any noise, absent from the 
asteroseismology measurements, we compute the rms of syn- 
thetic RV measurements for different time spans. To be consis- 
tent with the observational strategies we will use afterwards, we 
calculate the rms for 1 measurement per night of 15 minutes 
for 10, 50, 100 and 1000 consecutive nights. For each star, each 
measurement of 15 minutes is divided in several exposures ac- 
cording to the exposure time of the asteroseismology run (see 
Table[T]i. The difference between the synthetic RV rms for each 
time span and the rms of the real RV is less than 10% for all 
the stars. This coiTesponds to less than 25 cms"' which is sig- 
nificantly below the long term precision level of HARPS (60-80 
cms"', HARPS team, private communication). This proves that 
computation of synthetic measurements over long period does 
not introduce any bias. Therefore, such measurements contain 
only short term perturbations (oscillations, granulation phenom- 
ena, photon and instrumental noise and a little bit of activity) 
even if the total time span is much longer. 

Table [3] gives the rms for the real RVs, rmsdaia and for the 
synthetic RVs for 100 consecutive nights, nnSsy„jor- We also 
calculate the rms expected for each noise contribution as well 
as the high frequency noise derived using the power spectrum. 
This high frequency noise corresponds to the median level of the 
power spectrum for frequencies higher than p-modes. We notice 
a good agreement between the rms calculated using the fitted 



* Using the phase calculated by the periodogram is not what we want 
because it will reproduce the real data, setting the RVs at during the 
observation gaps. 
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Table 2. Parameters for the spectrum fit (SG - supergranulation, MG = mesogranulation, G = granulation). Asc, Amc^ Ac, and 
the constant are expressed in power density, ^^^^ 



liHz 



Star 


AiG 


Bsc [h] 


Csa 


Amg 


Bmg [h] 


C-MG 


Ag 


Bo [min\ 


Cg 


Al 


r [mHz] 


vo [mHz] 


constant 


ySHyi 


0.055 


24.3 


4.3 


0.021 


3.4 


4.0 


11.610-' 


72.8 


5.0 


18.5 10-3 


0.17 


1.0 


2.6 lO-* 


yu Ara 


0.029 


13.0 


6.0 


0.027 


3.4 


5.0 


1.1 10-3 


43.8 


4.5 


3.2 10-3 


0.26 


1.9 


5.1 10-* 


Q-CenA 


0.027 


7.4 


3.1 


0.003 


1.2 


3.9 


0.3 10-3 


17.9 


8.9 


2.6 10-3 


0.36 


2.4 


1.4 10-* 


rCeti 


0.027 


6.7 


2.6 


0.002 


1.2 


8.9 


0.3 10-3 


18.5 


19.8 


0.3 10-3 


0.75 


4.5 


1.7 10-* 


aCenB 


0.002 


12.0 


4.8 


0.001 


0.7 


4.4 


0.1 10-3 


8.9 


7.5 


0.2 10-3 


0.68 


3.9 


0.5 10-* 



Table 3. Comparison of the dispersion of real RVs, rmsdata, with the synthetic data dispersions for all type of noise, rmSsY„jot, and 
for each type of noise individually (SG = supergranulation, MG - mesogranulation, G - granulation, const = instrumental and 
photon noise). The dispersion of synthetic measurement is calculated on 100 consecutive days, with 1 measurements per night of 
15 minutes, rmshigh frequency corresponds to the noise calculated in the power spectrum at high frequency, after p-modes. 



Star 


rmSdata 




rmSjy„,sG 


rms jy„, MG 


rms,,„,,G 


rmS_yy„_ p-modes 


rms^^y,,^ const 


rmS/i/g/, frequency 




[m s-'] 


[m s-'] 


[ms-'] 


[m s-'] 


[m s-'] 


[ms-'] 


[ms-'] 


[ms-'] 


ySHyi 


2.60 


2.88 


0.55 


1.05 


1.26 


2.23 


0.94 




/J Ara 


1.97 


1.90 


0.56 


1.05 


0.46 


1.11 


1.00 


1.08 


aCen A 


2.03 


1.89 


0.77 


0.62 


0.40 


1.20 


0.99 


1.1 


rCeti 


1.47 


1.42 


0.85 


0.43 


0.34 


0.52 


0.87 


0.96 


aCenB 


0.9 


0.88 


0.15 


0.48 


0.32 


0.39 


0.52 


0.52 



constant and using the power at high frequency, which guaranty 
a correct estimation of the photon and instrumental noise. 



As already pointed out in llVLavor et al.l (l2003l) . we observe 
for dwarf stars a trend between the level of RV variation and 
spectral type. If we look at the main sequence stars in our sam- 
ple {a Cen A (G2 V), Ara (G3IV-V), r Ceti (G8 V) and a Cen B 
(Kl V)), we notice a decrease of the level of RV variation from 
G to K dwarfs. Figure |3] (top panel) illustrates well this trend. 
A similar behaviour can also be pointed out when comparing 
the RV variation with surface gravity, log g. In this case, all the 
stars follow the trend, even p Hyi (G2IV) which is a sub giant. 
Thus, it seems that the RV variation is dependent on the spectral 
type and the evolutionary stage of the star The expected results 
on averaging the stellar noise will then strongly depend on the 
characteristics of the star considered. 

Since the fit is the sum of each noise contribution, see Eq. |5] 
we can extract each noise independently from the VPSD and cal- 
culate the corresponding rms. Except for supergranulation, we 
can see in Table [3] that the noise level of granulation phenom- 
ena, as well as oscillations, decreases when we go towards late 
spectral type stars (the stars are ordered from early spectral type 
to late spectral type and from evolved star to non evolved). The 
calculation of the noise rms induced by supergranulation is very 
sensible to the parameters B and C as well as the minimum fre- 
quency of the spectrum, which explains why we do not see any 
trend between the noise level of supergranulation and the spec- 
tral type. The effect of noise amplitude in function of the spectral 
type and/or evolution is better seen in Table |2]by comparing the 
values for A^c, Amc, Ac and A/,. In this cas e, we can s ee th e 



Kieldsen & B edding (i 19951 



O' Toole et aL(.2008.) already 



dependency of superg r anulat ion . 
IChristensen-DalsgaardI (l2004 andC 
showed the dependency between amplitude of stellar noise and 
spectral type and evolution for oscillation modes. At our knowl- 
edge, it is the first time that this dependency is shown for gran- 
ulation phenomena. In conclusion, early K dwarfs have a total 
noise level lower than G dwarfs or sub-giants, making them bet- 




4 5 8 10 

spectral type (G0 = 0, G5 = 5, K0 = 10...) 




Fig. 3. Top panel: RV variation, rmsdata, in function of the spec- 
tral type. The fit is done on dwarf stars (circles) without taking 
into account the sub-giant (triangle). Down panel: RV variation 
in function of the surface gravity, log g. The linear fit is done on 
all the sample. 



ter candidates for the investigation of exoplanets using precise 
radial- velocity measurements. 

All the stars in our sample are very bright, below magni- 
tude 5.2 (see Table [TJ, which induce a very small photon noise. 
For yu Ara which is the fainter one, the photon noise estimated 
for a single exposure of the asteroseismology run is below 34 
cms-'. Comparing this value with the one found for the photon 
and instrumental noise (see imSsy„ const in Table |3]l, we see that 
the precision is limited by the instrument. We can separate our 
sample in 2. On one side aCenB which as a very low photon 
and instrumental noise, and on the other side, the other stars. 
Measurements for a Cen B are very recent and a lot of work as 
been done to reduce the effect of guiding noise, which explains 
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the rms difference between the 2 samples. In the sample with 
high rms, all the stars seem to have a similar instrumental noise, 
around 0.9 to 1 m s ' . Several factors can be responsible for the 
small rms difference. For example, the quality of the sky plays an 
important role in the instrumental noise, mainly guiding noise. If 
we look at a star with a very good seeing, smaller than the spec- 
trograph fiber diameter (1 arcesec), the point spread function of 
the star will not cover the entire fibre. This will introduce a noise 
when the star will change of position on the fiber because of im- 
perfect scrambling. This noise can be average for long exposure 
times, but not for high observational frequency as requested by 
asteroseismology. 

3.4. Correcting instrumental noise for long period 

As shown in last paragraph, our RV measurements seem instru- 
mental noise limited. This noise is coming mainly from the tele- 
scope guiding system which is not optimized to observe stars us- 
ing very short exposure times. For planet surveys, we normally 
use exposure times of 15 minutes which allow a good averaging 
of the guiding noise and should give smaller instrumental noise. 
Moreover, estimating the instrumental noise with asteroseismol- 
ogy measurements gives a good approximation of the noise on 
short time scales, but not on the long time scales required for 
planet surveys. For long time scales and 15 minutes exposure 
time per measurement, the instrumental noise on HARPS is es- 
timated to be between 0.6 and 0.8 ms"' (HARPS team, private 
communication). To take into account long term instrumental ef- 
fects and remove short term perturbations due to guiding, keep- 
ing the same contribution from oscillations and granulation phe- 
nomena, we readjusted the level of the fitted constant to 0.8 m s ' 
for all the stars and recalculated the synthetic RVs. This value of 
0.8 m s"' has been taken to be conservative. 



4. Considered observing strategies 

Using the synthetic RVs for each of the considered stars, the 
goal is now to estimate, for different observational strategies, the 
variation level of binned RVs induced by stellar noise. 

4.1. Present HARPS-GTO observational strategy 

In order to be useful in a practical way, simulations have to 
take into account a realistic access to telescope time. A suitable 
approach is provided by the scheduling of the hi gh-precision 
HARPS-GTO program. A star like HD 69830 (see iLovis et"aD 
[2006, Lovis et al. 2010 in preparation for additional points) as 
been followed using an ideal calendar of 1 measurement per 
night on 10 consecutive nights per month and over more than 
4 years. In the present study, we will use this ideal calendar on 4 
years and remove 4 month per year, since most of the stars dis- 
appear from the sky during this period. In addition, we suppress 
randomly 20 % of nights, to take into account bad weather or 
technical problems. Such a calendar represent thus 256 nights of 
measurements over a time span of 4 years. 

For each star in our sample, we calculate the RV expected 
for each date of the calendar using an exposure time of 15 min- 
utes, corresponding to the present HARPS-GTO observational 
strategy. This strategy will be referred hereafter as the IN strat- 
egy. Then we make bins of 2, 5 and 10 consecutive nights and 
calculate the rms of the binned RVs, rmsRvh, in function of the 
binning in days. 



Table 4. Characteristics of the different observational strategies 



Measure 


Number of 


Exposure 


Total time 


spacing [h\ 


meas ./night 


time [min} 


per night [min] 


2 


3 


10 


30 


2 


3 


20 


60 


5 


2 


15 


30 


5 


2 


30 


60 


24 


1 


15 


15 


24 


1 


30 


30 



The result for each star is illustrated by the heavy solid black 
lines with circle markers on Fig.|4] For example, in the case of 
aCenB, we can see that with a binning in days of 10, we can 
reduce the RV variation, mainly due to stellar noise, down to the 
level of 0.38 ms to be compared with 0.94m s~' in the case 
without binning. 

4.2. Setting new observational strategies 

We have presented in the previous section the result expected 
when using the current HARPS observational strategy for the 
high-precision subprogram (1 measurement per night of 15 min- 
utes). The goal of our work is to determine if we can reduce 
even more the noise contribution by sampling better the time 
scales of the various intrinsic stellar noise sources. The follow- 
ing additional strategies have been tested: 3 measurements per 
night with 2 hours of spacing and 2 measurements per night 5 
hours apart. For the first strategy, we tested for each exposure 10 
and 20 minutes and for the second one, 15 and 30 minutes. In 
addition a strategy with 1 measurement per night of 30 minutes 
has been tested. For long-period planet investigation, we more- 
over need to follow the star as long as possible along the year in 
order to obtain a good coverage in phase. Also, because of air- 
mass limitations, we consider that stars can only be observed at 
most 5 hours per night under good conditions. This explains the 
choice of the selected approaches. Table|4] presents the different 
strategies used in the simulations. 

The strategy with 15-minutes total exposure time per night 
corresponds to what is done in the high-precision HARPS- 
GTO subprogram. The corresponding results will be useful for 
comparison with results obtained with the new considered ap- 
proaches. The different exposure times considered will give us 
a good feeling on how to average stellar oscillation effects (p- 
modes that have a typical time scale of a few minutes for solar 
type stars), whereas the multiple exposures during the night will 
tell us about possible dumping, after binning, of granulation phe- 
nomena effects at lower frequencies. 

5. Results for the selected observational strategies 

In this section, we focus on the improvement brought by the new 
strategies. First, we calculate the level of RV noise, nnsRvh, as 
a function of the binning in days. We use the same calendar as 
in the previous section, which will allow us to compare the new 
strategies with the present one used on HARPS. 

5.1. An efficient and affordable observational strategy 

The evolution of the RV variation, nnsRvb, as a function of the 
binning in days is shown in Fig. |4] In these plots, each type of 
line corresponds to the same number of measurement per night. 
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Fig. 4. Rms of the binned RVs, rmsRvb, as a function of the binning in days. The different considered strategies are 1 measurement 
per night (the 2 heavy continuous lines), 2 measurements per night with a spacing of 5 hours (the 2 dashed lines) and 3 measurements 
per night with a spacing of 2 hours (the 2 thin continuous lines). Circles represent a total observing time of 15 minutes per night 
(1x15 minutes), squares represent a total observing time of 30 minutes per night (1x30 minutes, 2x15 minutes, and 3x10 minutes) 
and triangles represent a total observing time of 1 hour per night (2x30 minutes and 3x20 minutes). 



For each of them, we have 2 curves, the lower one is obtained 
with an exposure time twice longer than the upper one. The dif- 
ference between these 2 curves is small, so doubling the expo- 
sure time does not strongly improve the results, although it dou- 
bles the total measurement cost (neglecting overheads). Thus, it 
appears that the exposure time is not the only parameter that can 
average out stellar noise. The frequency of measurements plays 
also an important role. 

Taking only the shortest exposure times, we are left with 
3 strategies: 3 times 10-minute measurements per night with 2 



hours of spacing between them (hereafter, 3N strategy), 2 times 
15-minute measurements per night separated by 5 hours (here- 
after, 2N strategy), and 1 measurement per night with 15 min- 
utes exposure time (hereafter, IN strategy). Comparing the 1 
measurement-per-night of 30 minutes strategy with the 2N and 
3N strategies allows to clearly see that with a similar observation 
time per night (30 minutes), splitting the measurements over the 
night in 2 or 3 blocks improve significantly the averaging of the 
stellar noise. The best among the considered strategies is the 3N 
strategy. This strategy gives values for nnsRvb in average 30 to 
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40 % smaller than the present strategy used on HARPS for high 
precision. 

As already pointed before, the RV noise due to stellar in- 
trinsic perturbations is lower for K than for G dwarfs. Within 
a given spectral type, we can also compare this value between 
a main sequence and a more evolved star (a Cen A (G2 V) and 
fi Ara (G3IV-V) orySHyi (G2IV)). For more evolved stars, the 
RV variation seems higher, though we need to be careful since 
we just have one target, /3 Hyi (/i Ara is between a dwarf and a 
subgiant). 

5.2. Strategy effects for the different stellar noise sources 

The impact of exposure time and frequency of measurements on 
each kind of intrinsic stellar noise can be obtained by consid- 
ering separately each noise contribution in the power spectrum. 
Since we fit each type of noise, we can extract independently the 
3 types of granulation phenomena noises and the noise coming 
from oscillation modes. Then, using the technique to obtain syn- 
thetic RV measurements, we can create some data containing 
each type of noise, apply different strategies on them and thus 
see the impact of exposure time and frequency of measurements 
on each kind of stellar noise. 

Without surprise, p-modes, causing oscillations of the solar 
surface over time scales smaller than 15 minutes can be ade- 
quately averaged by longer exposure timefl For slightly lower 
frequency noises, it is useful to increase the total exposure time 
as much as possible (~ 30 minutes). This will average the effect 
of granulation whose typical time scale is less than ~25 minutes. 
For the mesogranulation, a further increase of the exposure time 
is poorly efficient whereas observing the star more often dur- 
ing the night leads to better results after binning. The 3N strat- 
egy is the one providing the best results for mesogranulation, 
since the 2 hours of spacing between the measurements over the 
night manages well to sample the perturbation time scale of the 
phenomenon. Finally, for supergranulation at even longer time 
scales, taking 2 measurements with a spacing of 5 hours gives 
the best results. 

6. Comparison with real observation 

Our simulation only takes into account oscillations, granulation 
phenomena, instrumental and photon noises. Lower frequency 
variations, due to active regions, are supposed to produce a 
higher level of noise, even for stars in their minimum of activity 
(iMeunier et al.ll2010h . To check if our simulation, not consider- 
ing active regions, corresponds to a realistic case, we have to 
compare the level of noise simulated with long term jitter ob- 
served on other stars. 

Our simulation gives us the noise expected for a given obser- 
vational strategy. In order to compare this simulated noise with 
real observations, we have to use an observational strategy simi- 
lar to the one nowadays used on HARPS. We thus choose the IN 
strategy. Using this strategy on a Cen A and a Cen B with no bin- 
ning, the dispersion taking into account stellar, instrumental and 
photon noise is equal to 1.29 and 0.94ms respectively. These 
values can be compared to the rms of the residuals obtained for 
detected planetary systems. Table|5]gives the rms of the residuals 
for host star presenting a spectral type similar to a Cen A (G2 V) 
and Q-CenB (KIV). The obtained values are in good agreement 
with observations, which indicate that our simulation of stellar 

' We can take 1 measurement per night of 15 minutes or 3 measure- 
ments per night of 5 minutes and do some binning, the results are similar 



and instrumental noises are realistic and that stars with noise 
dominated by high and moderate frequency variations exist. 

7. Detection limits in the Mass-Period diagram 

In order to derive detection limits in terms of planet mass and pe- 
riod accessible with the different studied strategies, we calculate 
the false alarm probability (FAP) of detection using bootstrap 
randomization (Endl et al. 2001; Efron & Tibshirani 1998). We 
first simulate a synthetic RV set containing stellar, instrumen- 
tal and photon noises as we did in the previous sections. The 
calendar used (see Sec. 14. Il l is the same, with 256 nights over 
a time span of 4 years. Then we caiTy out 1000 bootstrap ran- 
domization^ of these RVs and calculate the corresponding pe- 
riodograms. For each periodogram, we select the highest peak 
and construct a distribution of these 1000 highest peaks. The 1 % 
FAP corresponds to the power which is only reached 1 % of the 
time. The second step consist in adding a sinusoidal signal with a 
given period to the synthetic RVs. We calculate the periodogram 
of these new RVs and compare the height of the observed peak 
with the 1 % FAP for the given period. We then adjust the semi 
amplitude of the signal until the power of the peak is equal to the 
1 % FAP. The obtained semi amplitude corresponds to the detec- 
tion limit of a null eccentricity planet with a confidence level of 
99 %. To be conservative, we test 10 different phases and select 
the highest semi amplitude value. 

In practice, we simulate 100 RV sets for each strategy, and 
only the 10 "worst" cases were considered in the end. In Fig.|5] 
we can see the mass detection limits, with a confidence level of 
99 %, for the IN, 2N and 3N strategies using HARPS. We note 
on each graph a peak of the mass detection limits at one year. 
This can be explain by the 4 months removed each year in the 
calendar to take into account that stars disappear from the sky. 
Signals with 1 year of period are not well sampled, which leads 
to a local raise of the detection limit. 

Among the 3 strategies studied, we notice that the use of the 
2N and 3N strategies allows us to decrease the level of detection 
limits by approximatively 2 Me, compared to the high precision 
strategy normally used on HARPS (IN strategy). This improve- 
ment is not negligible, since it reduces mass detection limits by 
approximately 30 %. The difference between the 2N and the 3N 
strategy is small, but as shown before in Sec 15. II the 3N strategy 
seems the most efficient one to average out the considered stellar 
noise. 

In Fig|5] stars are ordered from top to bottom and left to right 
from early spectral type and evolved to late spectral type and not 
evolved. Following this order, we notice that the level of mass 
detection limits is going down. Thus, o-CenB, which is a Kl 
dwarf star, has the smallest mass detection limits in our sample. 

If we only look to dwarf stars, aCen A (G2V), rCeti (G8V) 
and a Cen B (KIV), we notice that the level of the mass detection 
limits is going down when we go towards late spectral types. 
For example, for 100 days of period, the mass detection limit for 
the 3N strategy will be 3.5 Me, 2.5 Me and 2 Me for the G2V 
star, the G8V and the KIV, respectively. Thus, late spectral type 
stars have lower mass detection limits, which can be explained 
by 2 phenomena. The first one, obvious, is that early K dwarfs 
will be less massive than G dwarfs. Therefore, for a given mass 
planet, the RV signal will have a higher variation for late spectral 
type star Nevertheless, this effect is not sufficient to explain this 

Following the use of this well sampled calendar, we choose to cal- 
culate the FAP using bootstrap rand omization rather than per mutation 
since it is more efficient in this case dEfron & Tibshiranilll998h 



10 Dumusque et al.: Stellar noise and planetary detection 

Table 5. Spectral type, activity level and rms of the residuals, Rmso_Ci for known planetary systems discovered using HARPS 



Star 



Spectral type \og(R'^^) Rmso-clfns '] 



Reference 



HD69830 KOV -4.97 0.81 Lovis et al. f2006') 

HD40307 K2.5V -4.99 0.85 Mayor et al. (2009b) 

HD47186 G5V -5.01 0.91 Bouchv et al. (2009) 

HD4308 G5V -4.93 1.3 Udry et al. (2006) 

HD181720 GIV -5.01 1.37 Santos et al. (2010) 
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Fig. 5. Mass period diagrams for our star sample. We can see on each graph, the 3 strategies studied. 1 measurement of 15 minutes 
per night (continuous thick line), 2 measurements per night of 15 minutes 5 hours apart (dashed hne) and 3 measurements per night 
of 10 minutes 2 hours apart (continuous thin Une). 



large range of mass detection limit^ and the second phenomena, 
which is the most important one, is that early K dwarfs have a 
lower level of stellar noise, compared to G dwarfs (see Table[3] 
and FigHJi. 



' aCenB as a masse of 0.90 Mq, which is slightly smaller than 
aCenA, 1.09 Mq 



If we now look, in the same spectral type range, to stars 
with different evolutionary states, we notice also an interest- 
ing behavior p Hyi is clearly a G2 subgiant, jj Ara is a G3 star 
in the transition be t ween the dwarf and the subgiant branch 
(ISoriano & Vauclahl 1201 Oh . whereas a Cen A is a G2 dwarf. 
Studying these 3 different cases, we notice that the level of mass 
detection limit is higher for evolved stars. For example, the sub- 
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Fig. 6. Mass period diagrams for aCenB. Thick and thin lines 
represent the IN and the 3N strategy, respectively. Plain lines 
corresponds to the simulated calendar and dashed lines to the 
real calendar ofHD 69830. 



giant p Hyi has mass detection limits 1 .5 times higher than the 
dwarf a Cen A. 

Following these results on dwarf and sub-giant stars, it seems 
clear that to find the lowest mass planets using RVs, we should 
look to early K dwarfs. Such stars present the lowest level of 
stellar noise and thus, the lowest mass detection limits. Besides 
this, early K dwarfs are also interesting because their habitable 
zones, around 200 days, are closer than for early G dwarfs. For 
stars like a Cen B (Kl V) (see Fig.|5]l, the simulations predict that 
we could find planets of 2.5 in the habitable zone, using the 
3N strategy. 

As already pointed before, due to a short time span of the as- 
teroseismology measurements, stellar activity noise is not fully 
included in our simulation. Thus, the present limits give the RV 
signal that can be reliably detected when activity noise is neg- 
ligeable. In the presence of real stellar activity noise, the detec- 
tion limits will raise up. This will be the point in a forthcoming 
paper. 

8. Detection limits using a real HARPS 
observational calendar 

In order to check if the calendar we use is not too idealistic com- 
pared to a real one (see section 14.11 ). we compute the detection 
limits using the real calendar of one of the most followed star 
using HARPS, HD 69830. This real calendar regr oups a total 
of 15 7 nights over an observing span of 1615 days dLovis et alj 
12006. Lovis et al. 2010, in preparation). Figure |6] shows the dif- 
ference in mass detection limits between the calendar used in 
this paper and the real one for HD 69830. At 100 days of pe- 
riod, the detection limit for the 2 calendar are 2 and 2.5 M®, re- 
spectively. This small increase of 25 % is only due to the total 
number of night present in each calendar. The semi amplitude 
in RV goes like the square root of the total number of measure- 
ment. Therefore, passing from 256 nights to 157 will increase 
the detection limit by 28 %, very close to what gives the simula- 
tion. Applying the 3N strategy on a K dwarf, such as HD 69830, 
should lead to detection of 3 M® planets in their habitable zone. 
Thus, by using a real HARPS observational calendar, we can 
reach similar detection limits. A real observational calendar de- 
pends not only on bad weather and instrumental problems, but 
also to the following time allocated for each star. The present 
strategy on HARPS is to follow a lot of targets, which reduce 



the number of measurement per star In searching for the small- 
est planets, just a very small sample of stars should be selected. 
This would allow us to increase the number of measurement 
compared to a real observational HARPS calendar. It is why the 
simulations, done in the previous section, use 256 nights of mea- 
surements on 4 years rather than 157. 

9. Concluding remarks 

The present high-precision HARPS observational strategy (1 
measurement of 15 minutes per night, IN strategy in the pa- 
per) manages well to reduce stellar oscillation noises, but is not 
optimized to average out granulation phenomena (granulation, 
mesogranulation and supergranulation). These latter sources of 
noise perturb radial velocity measurements on time scales up to 
1.5 days. Increasing the number of measurements per night al- 
lows us to reduce the effect of granulation phenomena, after bin- 
ning, and consequently improve the planetary detection limits. 
The best tested observational strategy, 3 measurements per night 
each 2 hours, with an exposure time of 10 minutes each (3N 
strategy in the paper), gives detection limits on average 30% 
better than the present HARPS high-resolution strategy. This im- 
provement is really due to the 3 measurements per night and not 
to the total observational time (which is doubled compared to the 
present HARPS observational strategy). Indeed, we have shown 
that 1 measurement per night of 30 minutes is not efficient. 

Our work also suggests a trend between the level of the stel- 
lar noise considered in this study and the spectral type for dwarf 
stars. Early K dwarfs have a total noise level lower than early 
G dwarfs. A trend between evolved and non evolved stars is 
also revealed, non evolved stars showing a lower level of stel- 
lar noise. Therefore, early K dwarfs seems to have the lowest 
level of stellar noise of our sample and consequently the lowest 
planet mass detection limits. In addition, the habitable zone is 
closest in early K than in early G dwarfs, making early K dwarfs 
the most promising targets to search for very low mass planets. 

Following the conclusion of the last 2 paragraphs, we applied 
the 3 measurements per night strategy on o-CenB (KIV). For 
the habitable zone of this planet, approximatively 200 days, our 
simulation shows that we would be able to detect planets of 3 
M® using an existing HARPS observational calendar Therefore, 
granulation phenomena and oscillation modes will not prevent 
us of finding Earth like planets in habitable regions. 

However, there are 2 different limitations to our present 
study. First of all, due to only short time span measurements, we 
do not take properly into account stellar activity noise sources, 
which can be important. The detection limits we obtain are only 
valid for stars without significant activity related phenomena, 
such as spots and plages. We have shown that such stars exist 
and therefore activity will not be a problem for some targets. A 
second limitation could come from the fact that we only simulate 
the presence of one planet per star The majority of small planets 
discovered nowadays are in multiple systems, which complicate 
the analysis. More data is simply needed in these cases to allow 
us to fit all the free parameters. 

Next generation spectrographs, such as ESPRESSO 
(http://espresso.astro.up.pt/) and CODEX (e.g IPasquini et al] 
2008), will reach better levels of precision and stability, most 
likely leading to a reduction of detection limits near 1 M®. 
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